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FOREWORD 


This  paper  was  prepared  at  the  Electronics  Sciences  Laboratory  as  part  of  the 
Independent  Research  program  of  LMSC.  It  is  the  second  of  a  series  of  contributions 
toward  finding  appropriate  diffraction  pattern  manipulation  techniques  to  improve  the 
discrimination  capability  of  aberrant  optical  systems  for  off-axis  objects  and  deals 
with  the  diffraction  in  astigmatic  systems.  Further  reports  will  be  concerned  with 
other  types  of  third-order  aberrations . 
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ABSTRACT 


On  the  basis  of  the  general  solution  of  the  diffraction  integral  for  a  rectangular 
aperture  in  the  presence  of  geometrical-optical  aberrations,  obtained  in  Part  I  of  this 
series,  this  paper  deals  with  the  effect  of  primary  astigmatism  upon  the  diffraction- 
pattern  intensity  distribution. 

The  diffraction  pattern  of  a  stationary  rectangular  aperture  is  described  by  means  of 
a  conveniently  chosen  Cartesian  reference  system.  We  consider  two  orthogonal 
aperture  orientations  and  show  that  the  intensity  distribution  depends  on  the  orientation 
of  the  aperture's  symmetry  axes  with  respect  to  the  meridional  plane  defined  by  the 
optical  axis  and  the  geometrical  image  point.  The  intensity  distributions  in  special 
cases  of  interest  are  indicated  for  each  of  the  two  aperture  orientations. 

In  the  case  of  a  slit  rotating  about  the  center  of  the  entrance-pupil  plane,  i.  e.  ,  in  the 
case  of  arbitrary  slit  orientation  with  respect  to  the  meridional  plane,  the  Cartesian 
reference  system  is  one  which  rotates  synchronously  with  the  aperture.  We  derive  a 
general  expression  for  the  diffraction-pattern  intensity  distribution  in  this  rotating 
coordinate  system,  and  obtain  its  mean  value  throughout  half  a  slit  revolution  by 
integration.  Again  we  indicate  special  cases  of  interest. 

Finally,  we  determine  the  overall  diffraction -pattern  intensity  distribution  of  the 
rotating  aperture  with  respect  to  a  space-fixed  coordinate  system.  The  accumulated 
intensity  at  any  arbitrary  point  of  a  space-fixed  observation  plane,  as  experienced 
for  example,  by  a  photographic  plate,  may  be  obtained  by  integration.  We  expect  that 
the  accumulated  intensity  at  the  Gaussian  image  point  due  to  the  continuously  rotating 
slit  has  a  greater  predominance  over  the  accumulated  intensity  at  a  point  in  the 
surrounding  field  than  has  the  intensity  at  the  geometrical  image  point  over  that  at  the 
field  point  when  the  rectangular  aperture  is  stationary. 
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C 

U(u, v  ,  w) 
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general  coordinates;  refer  to  Cartesian  reference  systems 

distances  of  the  geometrical  image  point  from  the  optical 
axis 

polar  coordinates  of  the  Gaussian  image  point  in  a  rotated 
reference  system 

dimensionless  coordinates  of  an  aperture  point  or  a  point  on 
the  Gaussian  reference  sphere  centered  at  the  geometrical 
image  point  and  passed  through  the  center  of  the  exit  pupil 

dimensionless  parameters  representing  width  and  length 
of  the  rectangular  aperture  respectively 

dimensionless  optical  coordinates  of  an  observation  point 
in  the  neighborhood  of  the  Gaussian  image  point 

I**1  power  of  aberration  function 

coefficient  in  power  series  representation  of  ,  n,  m 
being  running  indices,  £  a  constant  index 

Seidel’s  aberration  coefficient  of  primary  astigmatism 

value  of  the  disturbance  at  the  observation  point  ( u ,  v  ,  w ) 

light -wave  amplitude,  assumed  to  be  substantially  constant 
over  the  wave-front  portion  which  approximately  fills  the 
aperture. 

wavelength  of  monochromatic  light 
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Vj'm(u,v)  =  (-j)n"m  VJ"m(u,v) 
wJV.w)  =  (-j)mW^(u,w)  i 

Jr+ 1/2^* 

gJ(z) 


I(u,  V ,  w) 

yo.o.o) 

i(u  ,  v  ,  w) 


i(u  ,  v  ,  w) 
i(u.  vQ,  w 


o 


I( u  »  V  ,  W)/IQ 


1  * 

-  /  i(u,  V  ,  w)  d <p 
o 


<f> 

*0 


wave  number 

complex  solution  functions  occurring  in  the 
expanded  form  of  the  disturbance  U(u  ,  v  ,  w)  ; 
l  >  n  ,  m  being  running  indices 
th 

( r  +  1/2 )  order  Bessel  function  of  argument 
f  ,  r  being  an  integer 

solution  function  of  argument  z  ,  occurring  as 
factor  in  Bessel  function  series  expansion  of 
m  (or  W™  ),  r  ,  k  being  running  indices, 
defined  on  page  2-2 

constant  nominator  and  (common)  denominator 
of  coefficients  occurring  in  Bessel  function 
series  expansion  of  Gf(z  )  ,  derived  in  Refs.  1 
and  2,  p  being  running  index 
intensity  at  observation  point  ( u  ,  v  ,  w ) 

intensity  of  the  geometrical  image  point 
normalized  intensity  at  the  observation  point 

mean  intensity  at  relating  point  (u,v,w) 

instantaneous  normalized  intensity  at  space- 
fixed  observation  point  (u  *  VQ  »  WQ )  due  to 
rotating  aperture 

angle  of  initial  orientation  of  the  rectangular 
aperture  with  respect  to  the  meridional  plane 


<t> 


angle  of  the  slit  rotation  about  the  optical  axis 


radial  distance  of  observation  point  from  the 
u-axis 
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rotating  slit 
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Section  1 
INTRODUCTION 


In  a  previous  paper  (Ref.  1)  a  general  solution  of  Fresnel-Kirchhoff’s  diffraction 
integral  was  derived.  This  integral  pertains  to  a  rectangular  aperture  and  represents 
the  disturbance, due  to  diffraction, at  a  point  in  the  neighborhood  of  the  geometrical 
image  of  an  object  point  in  the  presence  of  geometrical  aberrations.  However,  the 
general  expressions  obtained  for  the  disturbance  and  for  the  (three-dimensional)  inten¬ 
sity  distribution  do  not  readily  yield  a  precise  insight  into  the  combined  effects  of 
diffraction  and  aberrations . 

For  such  an  insight,  it  is  helpful  to  consider  single  aberrations  affecting  the  structure 
of  the  diffraction  pattern.*  The  most  important  single  aberrations  are  contained  in 
Seidel’s  third-order  aberration  function.  This  paper  deals  with  the  effect  of  primary 
astigmatism . 


♦According  to  the  theory  of  geometrical  optical  aberrations,  the  describing  function, 
in  general,  is  given  by 


$  =  £  *<2k> 
k  =  2 

whare  is  a  polynomial  of  order  (2k)  in  four  significant  variables.  Every 

group  of  order  ( 2k  -  1 )  splits  into  a  finite  summation  of  specific  kinds  of  single 
aberrations,  and  each  aberration  has  a  different  effect  on  image  quality. 
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Section  2 

THE  DIFFRACTION  INTEGRAL  SOLUTION  IN  THE  CASE  OF 
ASTIGMATIC  ABERRATIONS 

If  the  meridional  plane  defined  by  the  optical  axis  and  the  principal  ray  of  the  image¬ 
forming  optical  system  is  chosen  as  the  YZ-plane  of  the  Cartesian  coordinate  system 
used  for  describing  the  diffration  pattern,  the  aberration  function  for  primary  astigma¬ 
tism  is  represented  in  a  very  simple  form. 

th 

In  the  notations  of  the  previous  paper  (Ref.  l),for  the  £  power  of  the  wave  aberration, 
we  have 


L 

n  =0 


n 

m=0 


„n-m  m 


n-m,m 


where  the  parameters  ^  m  generally  depend  on 

•  The  off-axis  position  of  the  geometrical  image  point 

•  The  linear  dimensions  of  the  aperture 

•  The  radius  of  the  Gaussian  reference  sphere 

•  The  lateral  magnification  of  the  pupil  planes 

•  The  aberration  coefficients  of  the  optical  system 

For  the  case  of  primary  astigmatism,  we  have  specifically 


$ 


(4)  = 


-  CY 


*2  2  2 
1  "o'* 
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Therefore,  the  coefficients  are  given  by 

K*  =0  for  (n-m)^0,m^2l 

n-m.m 


The  disturbance  at  the  abservation  point,  with  the  "optical  coordinates"  u ,  v  ,  w  , 
takes  the  form 


U(u, v,  w) 


-  *iu  ft 


4u  v 

o 


£  A  Ok)1 


£  =  0 


<2£V£°<U-V)W?<U,W) 


where 


vj(u,  v) 


=  f  £  ('if  (2s  +  1} 


s  +  1/2 


/u  2\ 

I  2  M  (;o 


s  =0 


(!-l 


I72~  G.<v»‘o> 


20 

Wpu,  w) 


(-  l/  f  J  (-  j/  <2t  +  x> 


Jt+l/2l 


t  =  0 


(H) 


G>“o> 


~  JN 

=  y  (-  d8-^ 

"  D 


a  =  0 


Ns,<r  J2s-2j+1/2(vmo) 

o  "  .1/2 

;  <v^> 


t  +  £  *t2£  T  /w  v 

„2t .  v  _  v  /  1  ,t  —  r  Nt,T  J2t+2£-2r+l/2(W,/o) 

1  0  L.  D,  (wv  )v 

T=0  t  O' 


•The  coefficients  ^  and  are  available  on  request  (Ref.  2) 
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The  intensity  distribution  is  given  by 


I(u,v,w)  =  I 


oo 

(Jk)1  K*o ,  2£  V® < u .  v)W^(u,w) 


where 


If  the  rectangular  aperture  is  rotated  by  ?r/2  ,  the  two  slit  parameters  and  uQ 
are  interchanged. 


By  performing  this  interchangement,  we  see  that  for  both  cases  the  disturbance  at  a 
specific  observation  point  (u ,  v  ,  w),  as  well  as  the  intensity  distribution, are  com¬ 
pletely  different. 

The  values  of  the  disturbance  function  and  the  intensity  distribution  in  three  charac¬ 
teristic  orthogonal  planes  are  believed  to  be  of  special  interest.  These  planes  are 

•  The  Gaussian  image  plane,  defined  by  u  =  0 

•  The  meridional  plane,  defined  by  v  =  0 

•  The  tangential  plane,  which  is  perpendicular  to  the  meridional  plane 
and  is  defined  by  w  =  0 


Since  we  have 


lim  Jr+l/2(f)  [2  6° 

£-0  ^172  \tT 
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where  6°  is  the  Kronecker  symbol,  we  obtain  for  u  =  0  nonzero  contributions  only 
from  terms  for  which  s  -  <r  =  t  =  0  .  For  the  meridional  plane  orientation  of  the 
slit's  longitudinal  axis,  this  yields, 


V°(0,v)  = 


G?) 


2 1 

W~(0,  w) 


N21  J, 


,  ffv  ,  ,,T  ‘■'o.T-a-2^1/2(WI'o> 

<■  1}  V2  2  1  "  ZTlUi 


T~0 


(WK0) 


i(0,v,w>  =  I(  0 ,  v ,  w )/ 1 


i{0,  v,  0)  = 


2  <-<*>  (CY1  -o)  ,'1>  D*«  (w„  1/! 

£-0  T=0  O  '  o' 


1=0  o 


i( o,  0,  w)  = 


i( 0,  0,  0)  = 


Y  1  ....£  (cy*2  2\l  f?  V  nT  No,T 

2  <? <)k>  (CYi  V  «  2  1  1(  Da 

1  =  0  r  =  0  o  o’ 


«  .  n2£ 

V  1  /  .ui/rv’2  2f  "o,i 

2  f7  (-  Jk)  (cv,  g 

jf  =  0  o 


If  we  replace  by  pQ  in  the  above  expressions ,  we  obtain  the  corresponding  results 
for  an  aperture  rotated  by  x/2  . 
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We  see  that  the  intensity  „t  the  center  of  the  dtffnmtton  pettem  I.  not  eqtud  to  muty 
an  ,  therefore,  differs  from  the  normalised  central  intensity  of  the  aberration-free 
Iffjaction  pattern  Furthermore,  from  interchmtgtag  po  ^  .  |, „,dM,  that 

Intensities  displayed  are  dependent  on  the  orientation  of  the  rectangular  aperture 
with  respect  to  the  meridional  plane. 


In  a  similar  manner,  for  the  meridional  plane  v  =  0  ,  one  obtains 


i(u,0,w)  = 


i(u,0,0)  = 


8-0  8  (%/2) 
2  A  <-  «  4  Vf  2  <-  »s  *  1 1 


00  /  o  V 

i/F 2  <-»‘ <»♦» 

t-0  B,a  /..  2,„11/2 


D‘  (-&) 

Finally,  for  the  tangential  plane,  w  =  0  ,  one  gets 


i(u,v,0) 


2  n  i-i)1  <*♦ i^wVK^  ' 

1=0  Dt  (u&f  1 


In  these  expressions,  the  corresponding  results  for  an  aperture  rotated  by  ir/2 
respect  to  the  meridional  plane  are  again  obtained  by  interchanging  ^  and  „ 


with 
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Section  3 

DIFFRACTION  BY  A  ROTATING  RECTANGULAR  APERTURE 

In  the  preceding  section,  we  have  noted  that  the  intensity  distribution  in  the  diffraction 
pattern  of  a  rectangular  aperture  in  the  presence  of  primary  astigmatism  depends  on 
the  orientation  of  the  slit  with  respect  to  the  meridional  plane  containing  the  geometrical 
image  point  and  the  optical  axis.  Of  course,  it  is  desirable  to  determine  the  shape  of 
the  diffraction  pattern  for  any  arbitrary  orientation  of  the  rectangular  aperture.  This 
determination  is  of  special  interest  when  we  consider  an  aperture  which  rotates  with 
constant  angular  speed  about  the  center  of  the  entrance -pupil  plane. 

The  solution  of  Fresnel-Kirchhoff's  diffraction  integral  is  usually  found  by  expressing 
analytically,  in  the  integrand,  the  optical  paths  involved.  The  intensity  distribution  is 
described  with  respect  to  a  Cartesian  reference  coordinate  system.  This  reference 
system  has  its  origin  at  the  (space-fixed)  geometrical  image  point  and  its  (positive) 
Z-axis  in  direction  of  the  light  propagation.  Although  the  orientations  of  the  X-  and 
Y-axes  may  be  chosen  arbitrarily,  it  is  convenient  to  place  them  respectively  in  the 
orthogonal  planes  defined  by  the  principal  ray  and  lines  drawn  parallel  to  the  aperture's 
lateral  and  longitudinal  symmetry  axes.  Without  restricting  the  generality,  we  may 
pass  the  slit's  symmetry  axes  through  the  center  of  the  entrance -pupil  plane. 

When  the  aperture  is  rotated  by  a  certain  angle  <f>  about  the  optical  axis,  the  diffraction 
pattern  rotates  synchronously.  However,  the  meridional  plane  and,  therefore,  the 
distance  of  the  geometrical  image  point  from  the  optical  axis  remain  unchanged.  For 
this  reason,  the  Gaussian  image-point  position  with  respect  to  the  reference  system 
is  expressed  in  polar  coordinates  as  follows: 

*  *  +  * 

=  Rj  sin  ip  ,  =  Rj  cos  <p 
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The  aberration  function  must  now  be  represented  in  its  general  form 


r  *  *  i2 

*  =  -  C  lXl  ^o  «  +  Y1  ■'o  *» 

-  l)1  (crJ2)  \(ho  sin  4>)  t  +  (uQ  cos  <*>)  i? 


and  its  Jt**1 


power  by 


i-  »'  «f  2  (?) 


(M0  8 in  <P) 


2i-\  .  .A  t2l-X  X 

(w  cos  0)  5  TJ 


X  =  0 


Hence 


■4-x.x  -  <-  •>' 


Kf(?) 


(po  sin  ^)2l~x  (|^0  cos  <*>)X 


This  yields 


o°  2£ 

U(u, v, w)  =  -  jeju  4po  i/Q)  £  j7  (jk)£  ^  ^-X.X^  X  (u.v)  (u,w) 


£=0  x=o 


where 


Vf-X(u,v) 


=  (-  J) 


2£  -  X 


12 


(-  j)“  (2s  +  1) 


Is+l/2l  2  V  „2t-J 
1>2  X  Gs 


s  =0 


(i“, 


o) 


(*«0) 
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«?(••»)  -  2 

t=0 


j)1  (2t  +  1) 


J 


(uv2  /2)1^2 


1/2 

G^lw-o) 


0a-x 

s 


('V 


28  +  2£  -  X  -  ( 1 ) 


(-  If 


2£  _  \ 

Ns  ,0-  J2B+2f-X-2<r+ 1/2^0  * 

,2£-X  .2  1/2 


(T~0 


S 


<V%>‘ 


GtX(Wl,o> 


2t  +  X  -  (1) 

2 

T=0 


J2t  +  X  -  2t+ 

(w.G) 


1/2^0  > 
1/2 


The  (normalized)  expression  for  the  intensity  distribution  takes  the  form 


i(u,v,w)  =  I(u,v,w)/lo 


*  21 

2f(-*<)  2  ■(?)«■. 


sin  <t>  )2<“*(  uQ  cos  <t>  )X  V^“X(u,  v)  W*(u,w) 


In  accordance  with  the  derivations  of  the  preceding  section,  for  <f>  =  0  ,  only  the  terms 

2£  2f 

containing  ,  and  for  <p  =  7r/2  ,  only  the  terms  containing  are  nonzero.  At 

every  observation  point  (u,v,  w ),  an  oscillatory  intensity  is  observed  as  <f>  varies 

between  0  and  2ir  . 

Note  that,  in  this  derivation,  the  specific  observation  point  chosen  does  not  remain 
space-fixed,  but  travels  along  a  circle  about  the  u-axis  in  step  with  the  aperture 
rotation. 
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Of  special  interest,  again,  are  the  intensity  distributions 

•  In  the  geometrical  image  plane,  u  -  0 

•  In  the  plane  v  =  0 

•  In  the  plane  w  =  0 

•  Along  the  u-axis ,  v  =  w  =  0 

•  Along  the  v-axis,  u  =  w  =  0 

•  Along  the  w-axis,  u  =  v  =  0 

•  At  the  diffraction  pattern  center  in  the  Gaussian  image  plane, 

u  =  v  =  w  =  0 


In  the  plane  u  =  0,  we  obtain  from 


2t-\-(l) 
2 


a -X  V  .  .  a  No,(T  Ja-X-2<7-H/2(V|‘o> 

''“o'  ■  Z  (_1)  Da-*  ^7 

(7=0  O  O 


x-(l) 

X  V  r  No.  r  JX-2t+1/2(WI'0> 
W  Z  <-  *>  D*  —  ,1 72 

T  =  0  O  O' 


cr  -  0  o 


him 

2  .A 


wVw,  =  |  (-  if  V  =  (-  7  wj(0,w) 

1  d  WI/J 

T  =  0  O  O 
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V^"A(0,v)  wj(0,w)  =  (-  l)1  vf ‘X(0,v)  wj(0,w) 

the  intensity  distribution 
i(0,  v,w) 


oo  21 

2  r.  (*<)'  2  (?)  ■*  t  ^  # 

1=0  £-0 


cos 


From  this  expression  we  obtain  the  distributions  along  the  w-  and  v-axes. 

2i  - 

that  in  both  cases  only  even  values  of  1  are  retained  in  the  functions  G 

X  ® 

Gq(0)  respectively.  Defining  these  even  values  by  X  =  2k  ,  one  has  a  = 

the  w-axis,  and  t  =  k  along  the  v-axis.  This  results  in 


21  -  2k 


of-V) 


(-  1) 


■  K  [2  ^  o  1  -  K 
r  d2£-  2k 


GoK<°> 


=  (-  1) 


v2k 

12  N,.. 

\ir  _.2k 


v‘f ~2k(o,o) 


2£-2k 

2e - 2k  .  -  JL - K  o,Jg-K 

i  d2£-2k 

O 


W2k(0,0)  =  (-  j)2K  (-  if  ^ 
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and  yields  the  intensity  distributions 


i(0, 0,w )  = 


2  r.  H2)‘  2  (2)  "o2 


2£-2k  2k  .2! -2k  ,  2k  . 
v  sin  <f>  cos  <f> 


1  =  0 


£  =0 


2£  -  2k 


x  (-  1)J 


i-K 


rv 

o,2-k  fif  X 
d2£-2k  v 2  2 


k  kt2k  t  .  . 

,  mt  No,r  J2k-2t+1/2(Wvo) 

<‘1)  - 7^72 - 

T  =  0  D0 


i(0,  v,  0)  = 


2  f  (*<f  2  &)  -o 

1=0  K  —  0 


2i-2K  2k  .  21- 2k  ^  2k  , 
Vq  sin  <p  cos 


k=  0 

i-K 


.  nr%‘ .  U^,T 

°*V\V  ’  -?-*•  ^ 


Hence,  at  the  geometrical  image  point, the  intensity  is  given  by 


i(  0, 0,  0) 


2  F.(->kcRf)‘  2  (L'K 

i=0  K=0 


^2£-2k  ^2k 

21 -2k  2k  _._2£-2k  . _ 2k  .  ^o.i-K  Ojk 

sin  4,008  4>—2 

O  0 


For  the  plane  v  =  0  ,  from 


of "2*(0) 


=  (- 1)' 


£  -  K 


^  m2*"2* 

4-^ 


.21 -2k 
8 
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8=0 


(*' 


0/ 


obtain  the  distribution 


i(u,  0,  w  )  = 


l~0  K=0 


2* -2k  2k  „^2f-2K 
0  o 


M„  C*  sin-  -  4  cob2k  <p 


X  W£K(u.w)  £  (-j)8(2s  +  l)-^±i^-S±iZ2_ 

Df  2*  /u  2\1/2 

8  U*/ 


s  -  0 


M 


^wvhere 
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/u  2 
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'  I 
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(-  1) 


,2k 


T  =  0 


(w^0) 


1/2 


S3  Similarly,  for  the  plane  w  =  0 ,  one  has 


i(u,v,  0)  = 


2  (”R 


1  =  0 


v  „2£-2k.  . 

x  v*  (u.v) 


2£-2k  2k  .  2{-2k  ,  2k 

sin  n  <f>  cos  4> 


k  =  0 


1-  j)  (2t+  1) 


t  =  0 


2k 

\t  +  K  Jt  4-1/2  1 
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M 
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where 


„2£-2  k.  . 
Vjj  (u.v) 


<-  it-*  i  f  (-i)*(»i*»»^4^  o?-* 

^  ^  /u  2\1/if 

U*w 


8  <V"o» 


s  =  0 


S+f-K 

I 

<r=  0 
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.  ,8-gNs%  J2s+2£ -2K-2g+l/2(v>io) 

3)  _2£  -  2k  “  ,1/2 


(v^o  > 


Along  the  u-axis,  finally,  the  intensity  distribution  is  determined  by 


i(u,0,0 ) 


00 

£  ^  .  . 

I  (“) 

8.  \  / 
1=0 

K=  0 

2£-2k  2k  .  21 -2k  ,  „  2ic 
li  v  sin  <p  cos  <p 


=o  „2£-2k  t  (±  „2\ 

F  S  t  ns  1 2s  +  n  a>a+trJL  8+1/z\2/°i 

V 2  Z  (  3).  (2S+1)  d2«-2«  /u  2\1/2 

(2  #*0) 


s  =0 


00  .,2k  t  /u  2\ 

[7  V  ltt  ,  Nt,t+jC  Jt+l/2U..t,oj 

V2  Z  (“3)  (2t+  X)  D2k  /u  2W2 

t=°  1  (2V 


As  indicated  on  page  3-3,  if  the  rectangular  aperture  is  rotated  about  the  optical  axis, 
a  specific  point  (u,  v,  w)  of  the  observation  plane  travels  along  a  circle  about  the 
u-axis.  Because  the  intensity  values  encountered  by  this  point  change  in  an  oscillatory 
manner,  the  determination  of  the  mean  intensity  value  which  the  traveling  observation 
point  experiences  along  its  circular  path  should  be  of  interest.  The  slit  orientation  is 
repeated  at  every  half  revolution;  consequently  we  have,  in  usual  notation, 

i(u,v,w)  =  —  f  i(u,v,w)  d<#> 
o 


3-8 


LOCKHEED  MISSILES  &  SPACE  COMPANY 


6-90-63-23 


For  convenience,  the  evaluation  of  this  integral  may  be  restricted  to  the  Gaussian 
image  plane,  since,  for  u  -  0  ,  the  intensity-distribution  function  contains  the 
imaginary  unit  only  as  a  factor  of  the  constant  |kCR*2j  .  Expressing  i(0,v,w)  In 
its  component  form,  one  obtains 


i(0 ,  v,  w)  = 


£  1>‘ xiy{kCRi2)J  2  (x) <*» '*  "o  <sta  <C°B  *>' 

£  =0  X=0 


,4£-X 


x  (0,v)  W,c(0,w) 


2 r 


+  j 


V  ,  ,.i  1  2\2J  +  1  V  /4f+2\  41+2-X  X,„,_iV4£+2- 

2  (_1)  (2£+TT!\kCRl  /  2  (  X  )Vo  Vo <8in^> 


4£+2 


£=  0 


X  =0 


*  va+i",<0,v*  Wa+i<°’w) 


In  performing  th<~  integration,  we  encounter  integrals  of  the  form 


r  ,  .  ,4(£  +  £' )  -  (X+  X' )  .  ^.X+X'  .. 

J  (sin  <p )  '  '  '  '  (cos  <p)  dtp 


and 


/  (sin  <M 


4(f  +  £’)  +  4-(X+X') 


( COS  0  ) 


X+X' 


d<£> 
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respectively.  These  integrals  are  known  to  become  zero  if  (X  +  X' )  is  odd.  For 
even  values  of  (X  +  X' )  ,  their  solutions  are  given  by 


ir 


(X  +  X'  )/2 

£  (-  if  ((Xt*',/2' 

p=o 


4(£  +  £')-(X+X')+2p  \  _ 1 _ 

2(£+£')+(X+X’  )/2  +  p/  24(£+£' )-  (X+X-  )+2p 


and 


(X+X'  )/2 

"  I 

p  =  o 


/  i  ,P  /(X+X')/2\  /  4(£  +  £')  +  4+(X+X')  +  2p  \  _ 1 _ 

\  p  /  \2(£  +  £')  +  2+(X+X' )/2  +  p/  24(£  +  £' )  +  4-(X+X')+2p 


respectively. 
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X=  0  X'  =  0 

(X  +  i')/2  . 
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2.  (  '  (  P  )  \2(l  +  V)+2-iK  +  \')/2*f)l 

P  *0 


,4«  +  l')  +  4-(X  +  X')  +  2p 
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This  formula  gives  a  rough  idea  about  the  mean  intensity  encountered  in  the  (rotating) 
diffraction  pattern  of  the  geometrical  image  plane  of  the  continuously  rotating  rectan¬ 
gular  aperture.  In  practice,  however,  this  diffraction  pattern  will  hardly  be  observable. 
Thus,  the  formula  derived  does  not  give  precise  insight  into  the  intensity  distribution  of 
the  overall  diffraction  pattern  of  the  rotating  slit. 

To  obtain  such  an  insight,  we  must  bear  in  mind  that,  as  the  slit  is  rotated  about  the 
center  of  the  entrance  pupil  plane,  every  space-fixed  point  (v0»w0)  °f  an  arbitrary 
observation  plane,  uq  =  constant  ,  successively  experiences  the  intensities  of  all 
points  along  a  circle  around  the  u-axis. 

For  convenience,  the  space-fixed  and  wQ-axes  may  be  placed  into  the  tangential 
and  meridional  planes,  respectively.  The  observation  point  may  then  be  described  by 
polar  coordinates 

v  =  p  sin  d>  and  w  =  cos  d> 
o  *o  ^o  o  *0  ^o 

If  we  now  rotate  the  slit  by  the  angle  <f>  -  and,  consequently,  the  reference  coordinate 
system  by  the  same  amount  -  we  may  determine  the  coordinates  of  the  (space-fixed) 
observation  point  by  the  transformation 

v  =  vq  cos  <p  +  wq  sin  <f>  =  pQ  sin  ( <pQ  +  <f>) 

w  =  -  vq  sin  <p  +  wq  cos  <p  =  pQ  cos  ( <t>Q  +  <p ) 

The  third  optical  coordinate,  u  =  uQ  ,  remains  unchanged.  These  optical  coordinates 
enter  the  intensity  distribution  function.  Of  course,  in  the  new  coordinate  system,  as 
in  the  preceding  derivation,  the  geometrical  image  point  is  described  by 

♦  *  *  * 

Xx  =  sin  <p  and  =  Rj  cos  <p 
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Thus,  we  obtain  the  instantaneous  intensity  distribution 


i(u,  v,  w) 


i  a  Kf  i  (?)  -r 

£  =  0  X  =  0 

x  (sin  (cos  X(U.P0  sin  (4>0  +  <t>))  w<  (u,po  008  + 


where 
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(2  w 
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Every  integrating  device  (as,  for  example,  a  photographic  plate)  accumulates  instan¬ 
taneous  intensities,  i.e. ,  throughout  half  a  slit  revolution,  it  measures  the  intensity 

i(u,v,w)tot  =  /  i(u, v,w)  d<*> 

0 

The  integration  involved  is  difficult  to  perform,  even  for  the  Gaussian  image  plane, 
u  =  0  .  We  expect,  however,  that  the  accumulated  intensity  at  the  geometrical  image 
point,  i(0,  0,  0)tot  ,  of  the  continuously  rotating  slit  has  a  larger  predominance  over 
the  accumulated  intensity  at  a  point  in  the  surrounding  field,  i(u,v,w)^  ,  than  has 
the  intensity  at  the  Gaussian  image  point,  i(  0, 0,  0 )  ,  over  that  at  the  field  point, 
i(u,v,w)  ,  when  the  rectangular  aperture  is  stationary.  Rather  than  perform  the 
integrations  for  a  variety  of  specific  field  points ,  we  prefer  to  rely  on  experiment  to 
reveal  the  expected  predominance .  Knowledge  of  this  predominance  will  indicate  one 
method  of  improving  the  discrimination  capability  of  an  optical  system  for  off-axis 
objects  by  appropriate  diffraction-pattern  manipulation. 
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